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Abstract 

We propose a g-deformed model of the anharmonic vibrations in diatomic molecules. We analyse 
the applicability of the model to the phenomenological Dunham's expansion by comparing with experi- 
mental data. Our methodology involves a global consistency analysis of the parameters that determine 
the g-deformed system, when compared with fitted vibrational parameters to 161 electronic states in 
diatomic molecules. We show how to include both the positive and the negative anharmonicities in a 
simple and systematic fashion. 

PACS: 33.15.Mt, 02.20.Uw, 31.15.Hz, 03.65.Fd, 02.20.-a 



1 Introduction 

Q-algebras were originally conceived as a means to solve the quantum Yang-Baxter equation but 
in the last decade a great number of applications have been found in diverse areas of physics, ranging 
from the deformation of conformal field theories |2j to optics j3| and nuclear and molecular spectroscopy 
On the other hand, algebraic models have been proposed and applied systematically in many 
fields, including in a relevant way the study of nuclei and molecules [Z1)|H|- This kind of approach has 
simplified the n-body problem dramatically and has given rise to numerous new insights, including, for 
example, a supersymmetric description of quartets of nuclei and the formulation of tractable models 
of polyatomic molecules [H]- In the latter case, these methods combine the use of Lie algebraic and 
discrete symmetry techniques, which describe the interactions and the global symmetry, respectively, of 
these systems. In particular, the su{2) algebra has been proposed as a basic algebraic structure, given 
its connection to the Morse potential |8j. This allows an algebraic treatment of anharmonicities, which 
become increasingly important as higher vibrational modes are excited, much in the same fashion as the 
harmonic oscillator algebra is used in connection with harmonic motion. 

While the Morse potential leads to a fairly good approximation to the spectroscopic properties of 
most diatomic molecules, it has some limitations. For example, its energy eigenvalues contain only 
quadratic corrections to the equally-spaced harmonic behavior, while additional, higher-order terms are 
usually required in the phenomenological Dunham expansion to produce an accurate fit to the observed 
vibrational energies. 

In itself, the Dunham expansion does not constitute an adequate theoretical framework, since no 
Hamiltonian, and hence no eigenfunctions, can be deduced from it. It is thus an interesting problem 
to see whether an algebraic framework can be constructed in such a way as to include these higher 
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order corrections, while maintaining its simplicity and providing a Hamiltonian and corresponding wave 
functions. 

In this paper we report on an investigation to extend the su{2) framework, by considering the q- 
deformation of the Morse Hamiltonian. We then test our ideas by comparing a series expansion of the 
g-Morse energy expression with experimental Dunham expansions for a large number of molecules. As a 
first step we restrict our analysis to test the consistency of this (/-extension with respect to the Dunham 
expansion. A subsequent and more delicate test should include the evaluation of the g-deformed Morse 
vibrational eigenfunctions and the calculation of dipole intensities in the same fashion as for the su{2) 
model. We wish to emphasise that, while other studies have considered the deformation of both vibrational 
and rotational degrees of freedom in molecules our approach is simpler and more basic, as we 

aim to generalise an established model for vibrational spectroscopy, that is not restricted to diatomic 
molecules but can be extended to polyatomic systems, by g-deforming the fundamental bond interaction. 
In addition, we test the g-deformation in a global way by comparing with the Dunham parameters of 161 
electronic states of diatomic molecules. 

We start the next section by discussing the traditional treatment of anharmonic vibrations via the 
Dunham expansion. 



2 Molecular Anharmonic Vibrations 

The phenomenological description of the vibrational energy of diatomic molecules in a given electronic 
state is given by the Dunham's expansion fTIT, traditionally written as: 

E{v) = hcUe (^"^ + ~ hcUJeXe (^^ + ^ + hoVeVe (^'" + ^ + " " " " 

where c is the speed of light in vacuum and v is the vibrational quantum number. The vibrational 
molecular constants We, WgXe and ojeUe are obtained by fitting the potential curve to the experimental 
spectral data, uJeXe <C Wg, uJeUe ^ WeXg. While the constant u>eXe is nearly always positive, the constant 
LOeUe can be positive or negative and is often very small ^]. The terms of quadratic and higher order 
with respect to {v + ^) in the expansion account for the anharmonic character of molecular vibrations 
and the constant Xg is called the anharmonicity constant. 

The expansion is part of the more general Dunham's description of the rotational-vibrational 
energy of molecules ^1 

E{v,j) = Y,m^(v + ^\ {j{j + i)r (2) 

l,m ^ ^ 

where J is the angular momentum and the coefficients yim are often called Dunham coefficients. Here, 
yio = hcLOe, 2/20 = —hcuJeXe, VsQ = hcuJeVe-, and the vibrational part is obtained by taking the terms with 
m = 0, i.e. by ignoring the rotational bands built on each vibrational bandhead. These parameters are 
essential since structural information is contained fundamentally in the vibrational spectra. Rotation- 
vibration interaction can usually be treated in a perturbative fashion. 

Dunham's expansion provides a convenient, empirical and model-free way of organising a large quan- 
tity of spectral data. It also provides a procedure for comparing data with the calculations arising from 
model potentials such as the Morse potential. This expansion, however, has the disadvantage that it does 
not arise from a Hamiltonian and hence does not provide wave functions for the vibrational states. 
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If the expansion Q is truncated to the quadratic term, one obtains essentiaUy the discrete spectrum 
of the Morse potential 



2^ 



The Morse potential [T^ . 



Em{v) = hcUJe ( + ^ ) ~ f W + ^ ) 1 • (3) 



VM{r-r,) = Djl-e-f^'-'^A^ (4) 



describes reasonably well the potential energy of diatomic molecules near the equilibrium position 
and accounts for the anharmonicity of the molecular vibrations. It is characterised by two parameters, 
De is the depth of the minimum of the curve and P (known as restitution constant in the harmonic 
approximation's models), /3 > 0. Using the Hamiltonian, 

Hm = -^V^ + Vm (5) 
Morse solved exactly the Schrodinger equation, and found the quantised energy levels 

Here is the reduced mass of the molecule. Comparing the coefficients in the expressions (jSJ and the 
well-known relations between the Morse parameters Dg and /3 and the Dunham coefficients are obtained 



De = ^ and /? = y^/xWeXg. (7) 

These relations are also known as consistency conditions jl4j . 

One significant feature of the Morse potential is that the number of bound energy levels is finite, i.e. 
V = 0,1,2, .. . [vm], where [vm] denotes the largest integer less or equal to vm- Indeed, the Morse energy 
^ has a maximum at vm, 

VM = U--A. (8) 



and thus the number of levels in the well is given by [vm] + 1- 

In what follows we will refer to the Dunham's expansion for vibrational energies truncated to the 
cubic term. It is convenient to write the energies in the form: 



l\ ^ iV ( \ 



E'(y) = E{v)/hj =\^ + -j-Xe\v + -j + ye [v + -J (9) 

where u) = 2iTCiJe- The truncation is natural as the experimental data show that each of the coefficients in 
the expansion is much smaller than the preceding one, and furthermore the contribution of terms beyond 
the cubic term is negligibly small for all known diatomic molecules and seldom used in practice 112j . 

The model we are proposing in the next section, aims at defining a Hamiltonian which leads to 
Dunham-like formulae for the molecular vibrations, using natural symmetry principles, starting from an 
appropriate deformation of the Morse Hamiltonian. 
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3 The Model and Its Approximations 



3.1 Hamiltonian 

We shall use a Hamiltonian that can accommodate in different approximations both the Morse energy 
^ and the Dunham expansion @. This Hamiltonian is given by 

H = a{J+J^ + J-J+) (10) 

where 7+ , J_ are the raising and lowering generators of the quantum group deformation of the algebra 
su{2), and a is a constant which we shall fix below. In order to be close to the usual formalism we shall 
realise these generators in terms of anharmonic g-bosons. Their algebra, known as quantum oscillator 
algebra HWg or q-boson algebra, has been introduced in ^1 ^1 and is a generalisation of the 
Heisenberg-Weyl algebra obtained by introducing a deformation parameter q. The algebra is defined by: 

aa^ — q^^a'^a = q^ , [n, a] = —a , [h, a^] = (11) 

where q is in general a complex number. This number is called deformation parameter since the boson 
commutation relations of the harmonic oscillator may be recovered for the value q = 1. The realization 
of Uq{su{2)) is taken from ^5j: 

J+ = a+[2j -n\g, = a, Jq = n - j (12) 
where the g-number is defined as: 

r 1 _ 9^ ~ 



q-q ^ 

In general j may be an arbitrary complex number, but in our context we shall take it to be real. These 
formulae indeed produce the standard relations: 

[Jo, J±] = ±J±, [J+, J-] = [2Jo]g (13) 

In general the basis of our system is determined by the application of the raising operator on the 
vacuum. The latter is denoted by |0) and is characterised by the standard properties - anninihilation 
by the lowering operator a and being eigenvector of the number operator: 

a|0) = , n|0) = i/|0) , (14) 

where may be an arbitrary complex number. The basis explicitly is: 

|n) = (at)«|0). (15) 

The action of the quantum group Uq{su{2)) on this basis is: 

Jo \n) = {n + v - j) \n) 
J- \n) = q^ [n\q \n — 1) 

J+|n) = [2j - V - n]q\n + I) . (16) 

For j a non-negative (half-) integer and = formulae (|12j) realise a unitary irreducible representation of 
Uq{su[2)) of dimension 2j -|- 1. To be close to this case below we shall suppose = 0. 
Substituting (|12j) in (|lUjl we obtain: 

H = a ([2j]q {[2]q [h]q + 1) - [2], [h]l q^'^ ) (17) 
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where e = ±1 and we have used 

aa^ = [h + l\q , a~^a = [n\q (18) 

and the (^-summation formula: 

[A + B]q=[A]qq^^ + [B]qq-^^ . (19) 
Using again ((T^ we rewrite (fTTj) to obtain: 

H = a ([n],[2i-n],[2], + [2i],) (20) 
The action of this Hamiltonian on our basis is: 

H \n) = £{n) \n) = a ([n], [2j - n], [2], + [2^],) |n) . (21) 



One motivation for H1U|) is that for q = I and choosing a = huj/Aj we get from ()21() (essentially) the 
Morse case: 

£{n)q=i = hw{n + i - n^/2j) (22) 
Note that this expression has a local maximum at n = j which is: 

(^.=l)max = nj),=l = ^^-^ (23) 

and the bound levels are below this value of j, i.e., we have the restriction: 

n = 0, 1, . . . , [j] (24) 
Comparing this with the Morse case we identify j with V]\j in Q which leads to the relation: 

— = 2j + 1 (25) 

This limiting case prompts us to choose in the general case: a = hu)/[2]q [2j]q and then we have: 

= (26) 

Then the eigenvalues are: 

c'f ^ _ c( M(n \ - Ng[2i-n]g , 1 _ sinh(n/p)sinh((2i-n)/p) ^ 1 

t W - t(nJ/(/^J - ^^.^^ +^2^^- ,i^i,(i/p),i^h(2j7p) +2cosh(l/p) ^''^ 

where we have introduced a parameter p so that g = e~^/^. 

This eigenvalue function is similar to the one given in formula (25) of though the approach in 
jOj is not the same to the present one, (a pair of g-bosons is used, the Hamiltonian is different), and 
furthermore, we note that in [Hj formula (25) does not follow (for q ^ 1) from the Hamiltonian given in 
formula (24). 

Note that in this form it is transparent that the eigenvalue function has only one extremum as the 
extrema condition is: 

sinh((2i - 2n)/p) ^ ^ 
psinh(l/p) sinh(2j'/p) 

and this is a local maximum at n = j (as in the undeformed Morse case (|22j) ). the maximal value being: 
which is, of course, a deformation of (|23|) . 
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3.2 The case of positive molecular constant ?/e 

One way to use this Hamiltonian is to suppose that q is real and close to 1. Then we expand to second 
order in 1/p : q = = \ — Ifp -\- l/2p^ + • • • and: 

q^ = l-m/p + m^/2p^ + . . . , [m]„ = ^™^("^/^) = rn{l + (m^ - l)/Qp^) + ■■■ (30) 

sinh(l/pj 

Further, one would consider the eigenvalues of h which are significantly smaller than p and apply (|30() 
for m ^ h. However, using this directly would not give exactly the Dunham expansion. In order to be 
closer to the Dunham expansion we expand the eigenvalue function (|27|) as a function of n around —1/2 
(as in 6 ) and obtain: 

£'{n) = (1 - sinh((2j + l/2)/p))/sinh(2j/7p))/2cosh(l/2p) + f(n + l/2)sinh((2j + l)/p)/p- (31) 



-(n + 1/2)2 (,osh((2j + l)/p)/p^ + 2(n + 1/2)^ sinh((2j + l)/p)/2>p^ j/sinh(l/p) sin\i{2j/p) 

This expansion is also appropriate in the region when n + 1/2 is much smaller than p. Note that we 
make no restriction on the values of j, an assumption which turns out to be justified. We also see that 
the coefficient of (n + 1/2)^ is positive which means that this expression corresponds to the case when 
the molecular constant (Dunham coefficient) is positive. The cases of negative molecular constant Ue 
will be considered below. 

We consider the expression in ()31() truncated to a cubic polynomial in (n + 1/2). The quadratic 
equation for its extrema has the following solutions: 



n± = I ( coth((2i + l)/p) ± V Mh((2i + l)/p))f - 2 ) - 1/2 (32) 



One is inclined to require that the discriminant be strictly positive: 

( tanh — ) 
V P J 



2j + l\^ 1 
2 



or 

2 j + 1 1 , , 

— < -ln(3 + 2V2 « 0.88 33 

p 2 

In that case there are two extrema. The extremum at n_ is a local maximum. Obviously the value n_ 
is a deformation of the Morse value j, and if j « p one can expand: 



(2i + 1)' 
12p2 



n-=j+ ' +••• (34) 



The extremum at n_|_ is a local minimum. The value n+ grows rapidly with p: 



p2 , 1 

n+ = — — - - 1/2 + const. ^ + 
2j + 1 



Thus, the eigenvalues n cannot be near and one cannot use the potential well around it. 

Thus, when the restriction H33|) holds, and similar to the Morse case, we shall be interested in the 
region: 

n = 0, l,...,[n_] (35) 
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However, j and p will be determined by the experimental data and we shall see that there are cases 
when (jSnj does not hold. In that case, the potential does not have a local maximum, i.e., there are cases 
when the Dunham expansion has a different behaviour than the Morse expansion. Yet we are interested 
in small values of n when the two potentials do not differ significantly. Thus, in these cases we shall be 
interested in the Morse region (|24|) : n = 0, 1, . . . , [j], though j shall be determined from the Dunham 
expansion. 

We would like now to check how the expansion (|3H) corresponds to the Dunham expansion, which 
can be done independently of whether ()33() holds or not. 

First, we take the ratio of the coefficients of the linear to cubic terms in (n+1/2) both in the Dunham 
expansion Q and in our 1)311 ). and we find: 

This means that p is determined from the value of i/e : 



Next we take the ratio of the coefficients of the linear to quadratic terms in (n+1/2) both in the Dunham 
expansion @ and in our (j^), and now we get: 

— = p tanh — (38) 

Xe p 

This is a deformation of the usual relation of the Morse model H25|) . and if j « p one can expand: 



1 (27 + 1)^ 

- = 2j + 1 - ^ % / 

Xe 3p^ 



Naturally, when p — > oo (t/e — ^ 0), the Morse case applies. 

More important is that (|38|) gives a test for the applicability of our model. Indeed, recall that the 
function tanh is restricted: tanhz < 1 for any real z. Thus, we have: 



1 > tanh = ^— = 1 

p pxe V 2,x^ 

Thus, we have the following restriction on our model from the experimental values Xg and ye'- 



Note that the limit |^ 1 with values below 1 corresponds to the limit j — > oo. Thus the breaking 
point ^ = I means j = oo. 

If the bound goes the other way around, i.e., % > |, then our model is not applicable (even if we 
consider complex j). 

In the cases when ()39p holds from H38|) we determine the value of j using the value of p from ()37|) , 

i.e., 



i = i U|;-tanh U|| -1 (40) 
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It is useful to write down the formula for n_ in terms of X(,,ye using (|37jl and (|38|): 



3ye \ V "^"^ 
from which it is clear that the restriction H33|) translates into: 



^- = ^ -1/2 (41) 



^ < i 

3?? 3 



2 < o (42) 



which is indeed stronger than (|39|) . but is not a restriction on the applicability of our model. 
In the case when ye/x^ << 1 we can use the expansion: 



and obtain from ()4U() : 



arctanh z = z + + 



3.3 The case of negative molecular constant ye 

In order to accommodate the situation when the molecular constant is negative we have to consider 
the deformation parameter g to be a phase (though not a root of unity), or equivalently to make the 
replacement: p ip, then q g = e'/^. Almost everything may be obtained from what we have in 
the case of real q by this replacement. In particular, one needs: 



sm{z/p) 
sin(l/p) 



sinhz/p —isin /p , [^]q ^Z7T7Z\ ^ cosh z/p ^ cos z/p , 



(l/p) sinha/p ^ —{1/p) sina/p, (l/p^) sinha/p ^ +{l/p'^) sina/p, 
sinh(l/p) sinh(a/p) i— > — sin(l/p) s'm{a/p) 
We start with the analogues of ((^ : 

H'^ = [n],[2j-h]J[2j],+l/[2]g, (44) 

and (EH) 

S'jn) ^^ ^in(n/p)sin((2j-n)/p) ^ 1 ^^^^ 

sin(l/j)) sin(2j7p) 2cos(l/j») 

This function is very different from (|27j) . in particular, it has an infinite number of extrema since the 
extremum condition is: 

sin((2j - 2n)/p) ^ ^ ^^^^ 
psin(l/p) sin(2j/^) 

i.e., the extrema are at: n = j + kpTr/2, k = 0, ±1, ±2, .... However, again we would be interested in the 
region restricted by the first positive local maximum: 

n = 0,l,...,[j] 

This is also consistent with our aim of fitting the Dunham expansion. The maximal value is: 

^ / sin(j7p) \ 1 

Vsm(l/p) J 2cos(l/p) ^ ^ 
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which is a deformation of 

The expansion of n around —1/2 is: 

£'_{n) = (l + sin((2i + l/2)/p)/sin(2j7p))/2cos(l/2p) + f (n + 1/2) sin((2j + (48) 



-(n + 1/2)2 cos((2i ^ i)/p)/p2 _ 2(n + 1/2)3 sin((2j ^ i)/p)/3p3 j/sin(l/p) sin(2j7p) 

i.e., this is suitable for the Dunham expansion with negative ye- We consider the expansion in (|48() 
truncated to a cubic polynomial in (n + 1/2). As such it has two extrema at the points: 



n± = I [^^{coi{{2j + l)/p)Y + 2 - cot((2i + l)/p)J - i (49) 

The value is negative and thus inaccessible. The extremum at n_ is a local maximum. 

Further, we take the ratio of the coefficients of the linear/cubic terms in (n+1/2) both in Dunham 
and in our (|15|l and we get: 

- = -^P^ < (50) 

Ve 2 

This means that p is determined from the value of ye '■ 



^"V 3^' ^^^^ 

Next we take the ratio of the coefficients of the linear /quadratic terms in (n+1/2) both in Dunham (jUj) 
and in our (|48l) and we get: 

— = p tan (52) 

Xe p 

From this we should determine the value of j using the value of p from (|5H) . however, further we have to 
distinguish whether > (which holds in most cases) or Xe < 0. 

• Xe > ^ tan ^ > 

Then we have for the value n_ : 



%e| 

"^lyel \ V Xe 



^- = I + ' ye<0, Xe>0. (53) 



It is a deformation of the Morse value, and if j << p we can expand: 



(2i + If 
I2p^ 



Thus, similarly to the Morse case we shall be interested in the region: 

n = 0, l,...,[n_] (55) 

The value of j is obtained from (|52() : 



9 



Since the function arctan is multivalued, in the last formula we take the value which is closest to the Morse 
value from H25|) . Analogously, we use this for the following expansion which is valid when —ye << Xg : 



• Xe < ^ tan ^ < 

In this case the relevant parameter is j' which is complementary to j w.r.t. pTr/2. Then, 
tan(2j + l)/p = tan (vr - {2j' + l)/p) = - tan(2/ + l)/p 
and instead of we shall use: 

J^=pta.^ (58) 

\Xe\ P 

Further, for the value at which there is maximum we have: 

n_ = I (^^/{cot{{2f + + 2 + cot((2j' + l)/p)^ - \ = 

''^ ' ^ 1 + ^ + (59) 



3|ye| V V X. 



e 



It is not a deformation of the Morse value, since it grows with p and will not be useful for our purposes. 
Thus, similar to a case above we shall be interested in the region: 

n = 0, 1, . . . , [j'] (60) 

The value of j' is obtained from (|58() : 



3 =\ \ \ arctan \ ] -l] , ye < , Xg < . (61) 




In the last formula we take the value which is closest to the Morse value w.r.t. j', i.e., j' = — ^• 
Analogously, we use this for the following expansion which is valid when —y^ « : 



4 Analysis of Experimental Data 

In the Tables below, we have calculated the independent parameters of the model p and j, using the values 
of the molecular constants published in jjjj for 161 electronic states of diatomic molecules. The molecules 
are listed in alphabetical order in the first column of the tables, the corresponding electronic states are 
given in the second column. For convenience, in the third and fourth column we display the values of the 
anharmonic molecular constants jijjf. and De = uJeVe/^e Calculated from the pubhshed data in 

|llj . The values of the parameters of the model p and j are given in the following two columns. In the 
last column, where appropriate, we give the value of 

In Table 1 we display the cases when the molecular constant Ue is positive and the values of the 
parameters p and j are calculated from (|H7j) and (jlQ]) respectively. In the cases when ye/ {xef' < 1/3, the 
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potential curve for a given electronic state has a local maximum and the maximal number of vibrational 
levels is determined by the value of [n_] calculated from H41|) . 

When 1/3 < ye/ (a^e)^ < 2/3, Dunham's potential curve truncated to the cubic term does not have a 
local extremum, which is indicated by "x" in the column for In these cases the maximal number of 
vibrational levels in the model is determined by the value [j] . 

As discussed in the previous section, when ye/ {xe)^ > 2/3, the model cannot be applied in its present 
form. The experimental data in show that this condition is satisfied for 30 electronic states of the 
following diatomic molecules: ^^As2 state X ^T,+, ^^B'^^Br state A ^U, 209Bii9F state A, ^^Cs state 

and X ^Hs state a ^T.j, Hg^^Cl state X 2^+, (Hg^H)"^ state X ^2+, (Hg^H)^ state X ^^Hn^R 
state X 39K2 state B ^Uu, ^Li^H state X (^^^Mg^R state A ^n^, ^^Mn^^O state A, ^SNaK 

state C ^n, i^N^^O state B ^Rr, ^^^2 state A 1^+, Pba state A, ^^Sea state (X) ^SSi^^^O state 

^ In, 28SiTe state E, YbCl state B {^R) and ^OZri^O state a (^S). 

In Table 2 we show the cases when the molecular constant ye is negative and the values of the 
parameters p and j are calculated from (|51() and (|56() respectively. Here, the model works without 
any restrictions. The number of vibrational levels is determined by [n_] calculated from (|53)) . The 
experimental data shows that when n_ ~ j, the value of j is at least three times smaller than the value 
of p. 

There are six electronic states A of the diatomic molecules ^^'^Cs^H, ^^^^K^H, ^Li^H, ^Li^H, 
^^Na^H and Rb^H in |llj . for which the values of both molecular constants Xe and ye are negative. The 
model can accommodate these cases and the values of the parameters p and j' , calculated from ()51() and 
(|61j) . are given in Table 3. The number of vibrational levels is determined by [j']. As we have pointed 
out the case < is not a deformation of the Morse potential. Thus, it does not behave like (deformed) 
oscillator models do. In spite of this, our model can handle it. 



5 Conclusions and Outlook 

In this paper we have considered a g-deformation of a general Hamiltonian constructed from the raising 
and lowering generators of the quantum group deformation of the algebra su{2) and have obtained an 
expression for the eigenvalues of this Hamiltonian which can be interpreted as the phenomenological 
Dunham's expansion of the vibrational energies of diatomic molecules. The expansion is truncated to 
the cubic term mainly for practical reasons as the experimental results suggest that the contribution 
of terms of higher order is often negligibly small and seldom used. We have formulated a model of 
the anharmonic vibrations of diatomic molecules which in different approximations leads to Morse or 
Dunham's results. The parameters of the model are obtained in terms of the well-known experimental 
molecular constants and this gives a clear test for the applicability of the model and its restrictions. 
The model can accommodate both the positive and the negative values of the anharmonic constant ye- 
Using all data available in we have tested the model in a global fashion for a large set of diatomic 
molecules and conclude that it fits well with the experimental data for all states, except for about 30 
electronic states for which the values of Xe are not much larger than the values of ye, as required by the 
model. In these cases the model could be expanded by adding to the Hamiltonian terms corresponding 
to the rotational energies in order to obtain a (^-deformed version corresponding to the general Dunham 
expansion @. Next, in order to test these ideas further, intensity data should be studied. This would 
involve a procedure where a consistent g-deformed dipole operator is defined and the wave functions for 
each molecule evaluated within our approximation. Finally, a more ambitious project would be to extend 
these studies to polyatomic molecules, following, for example, the methodology of reference [Hj. 
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Table 1: Parameters of the model for positive molecular constant ye 



molecule 


state 




Ue 


V 


j 


[n_] 

L J 


^^Al^H 


X is+ 


1.73217 X 10-^ 


1.545 26 X lO-'' 


65.68 


44. 53 


X 






1. 249 06 X 10"^ 


8. 086 14 X 10-^ 


90.80 


62.31 


X 


^Be^^O 


A 


7. 353 8 X 10~^ 


2.96180 X 10-^ 


150. 03 


112. 54 


X 




7.953 69 X 10-3 


1.502 70 X 10-^ 


210.63 


71.99 


81 






9.33829 X 10~^ 


3.965 3 X 10-*^ 


410. 03 


54.31 


55 






8 849 85 X 10~'^ 


3 886 73 X 10-^ 


130. 97 


84.90 


X 






3 1 78 95 X 1 0^-^ 


3 288 57 X 1 0-'' 


450 25 


194 16 


272 




X is+ 


6. 17701 X 10"^ 


1.413 36 X lO-'^ 


217. 18 


104. 02 


X 






8 655 24 X 10~^ 


8 38636 X 10"'' 


281.95 


60.87 


63 






2 7Q4 34 X 1 


3 710 61 X 10-^ 


42 3Q 


25 69 








2 688 73 X 10~2 


3 452 38 X 1 0~^ 


43 94 


26. 82 






X is + 


2.684 52 X 10"^ 


6. 59811 X 10-^ 


100. 52 


19.06 


19 






2 420 92 X 10~2 


2 822 Q5 X ^0~^ 


48 60 


30 02 


V 
A 




e'^S + 


2.71367 X 10~^ 


2.73991 X IQ-"' 


49.33 


23.33 


X 


a '^E + 


2.329 28 X 10"^ 


2.59916 X 10-'^ 


50.64 


31. 10 


X 




2.487 71 X 10^2 


3.816 78 X lO-'^ 


41.79 


40.66 


X 






1.962 79 X 10"^ 


1.43009 X 10-"^ 


68.28 


32.42 


X 


2tt 
H2 




2.21695 X 10-2 


1.843 71 X 10-^ 


60.13 


28.76 


X 




1.906 84 X 10-2 


1.802 91 X lO-'' 


60.81 


39.11 


X 






1.855 36 X 10-2 


3.606 96 X 10-5 


135.95 


28.01 


28 


1tt3tt 




2.242 67 X 10-2 


2. 369 73 X 10-"^ 


53.04 


31.95 


X 




2.19751 X 10-2 


2.305 91 X 10-^ 


53.77 


32.92 


X 






1.61321 X 10-2 


1. 158 80 X lO-'' 


75.85 


43.06 


X 


lH35ci 




1.740 95 X 10-2 


1.873 07 X 10-5 


188.66 


29. 16 


29 


(1h35ci) + 




2.46483 X 10-2 


2. 129 79 X lO-'' 


55.95 


25. 19 


X 


ijjiyF 




2.1764 X 10-2 


2.368 X 10-4 


53.06 


34. 43 


X 


Hg^H 




5.01733 X 10-2 


1.11842 X 10-3 


24.41 


13.49 


X 






2.272 73 X 10-2 


1.81818 X lO-'' 


60.55 


27. 40 


X 






1.46104 X 10-2 


1.623 38 X 10-5 


202.65 


30.26 


35 


KBr 


X is+ 


3. 030 3 X 10-3 


4. 761 9 X 10-*^ 


374. 12 


174. 08 


X 


KCl 




3.21429 X 10-3 


3.92857 X 10-*^ 


411.94 


164. 61 


X 






3. 301 89 X 10-3 


4.71698 X 10-^ 


375. 94 


208. 96 


X 


■'Li2 




6. 16168 X 10-3 


7.046 39 X 10-« 


307. 59 


89.76 


96 


^Li^H 




1.650 48 X 10-2 


1.16174 X 10-"^ 


75.75 


41.09 


X 






3.33333 X 10-3 


3.77778 X 10-*^ 


420. 08 


187.61 


X 


^%2 




6.12644 X 10-3 


3. 182 73 X 10-*^ 


457. 67 


84. 86 


87 


^'^NaBr 




3. 650 79 X 10-3 


2.539 68 X 10-*^ 


512.35 


152.33 


165 


^^Na^H 




1.682 31 X 10-2 


1.36495 X 10-"^ 


69.89 


43.46 


X 


^%al^^I 




2.622 38 X 10-3 


3.496 5 X 10-*^ 


436. 65 


293. 55 


X 


^«02 




7.63939 X 10-3 


3.45491 X 10-5 


138.91 


121.64 


X 


Rb^H 




1.51051 X 10-2 


8.00623 X 10-5 


91.25 


41.44 


X 






1.623 86 X 10-2 


1.136 95 X 10-4 


76.57 


42.01 


X 






4.868 64 X 10-3 


2.648 89 X 10-*^ 


501.68 


108.59 


112 


Sr^^F 




4.4673 X 10-3 


3.953 4 X 10-*^ 


410.65 


125.03 


136 


Zu^II 




0.0 343 


2. -175 7 X 10-^ 


51.89 


15. 99 


17 



12 



Table 2: ] 


Parameters of the moc 


el for negative molecular constant ye 


molecule 


state 






P 


i 


[n_] 


109Ag8lBj. 


B (^n-) 


2.4613 X 10"^ 


3.318 6 X 10-4 


44. 82 


16.00 


14 


luyAga.'^ci 


B (^n-) 


2.135 2 X 10"^ 


3.380 8 X 10-'' 


44.41 


17.53 


16 




^ (•■^n+) 


1. 286 X IC)-- 


3.908 8 X 10-3 


13. 00 


8. 07 


7 




A 


2.1534 X 10-^ 


1.773 2 X 10-3 


19.39 


10.90 


9 


2YA135C1 




9.712 X 10"^ 


4.800 4 X 10-* 


37.27 


22.30 


19 


^7ai19f 




1.032 9 X 10"^ 


2.272 5 X 10-** 


54. 16 


28.23 


25 






3.29789 X 10-^ 
1.870 7 X 10-2 


2.35393 X 10-3 
1.908 9 X 10-^ 


16.83 

186.88 


8.45 

25.53 


7 

25 






1.324 5 X 10"^ 


1.7615 X IQ-'' 


194. 54 


35.51 


34 


11B35C1 




1.340 7 X 10"^ 


1.179 2 X 10-** 


75. 19 


28.88 


26 






1.9064 X 10-2 


2. 395 X 10-4 


52.76 


20. 14 


18 




1.7245 X 10-2 


2.4288 X 10-4 


52.39 


21.40 


19 


(9Be^H) + 




1.002 6 X 10-2 
.01791 


2.5744 X 10-^ 
9.452 2 X 10-6 


160.92 
265.58 


44. 14 
27.02 


42 
26 






7. 743 5 X 10-'^ 
1.3261 X 10-2 


1.459 3 X 10-4 
3. 641 6 X lO-"' 


67. 59 
135.30 


36. 29 
33. 90 


32 
32 






5.650 3 X 10-^ 


1.969 6 X 10-^ 


1839.8 


87.72 


87 






2.275 92 X 10--^ 


3.585 21 X 10-6 


431.22 


170.86 


159 


XiE+ 


1.868 45 X 10-3 


1.343 29 X 10-5 


222.78 


130.53 


117 






3.92907 X 10-3 


7.57854 X 10-4 


29.66 


21.07 


18 






8.739 2 X 10-3 


2. 222 1 X 10-3 


17.32 


11.80 


10 




9.193 5 X 10-3 


2.833 5 X 10-4 


48.51 


27. 42 


24 




{A ^I^ 


9.295 8 X 10-3 


2.2715 X 10-5 


171.32 


47. 53 


45 


(x^ri)'^ 


7.312 3 X 10-3 


6. 728 X 10-5 


99.54 


46.36 


42 






6.688 2 X 10-3 


1.206 7 X 10-3 


23.51 


16.13 


14 


12^160 




6.700 7 X 10-3 


9.8876 X 10-5 


82.112 


43.34 


39 


— /10 ,1^; „\ — i — 

^i2QibQj + 


X ^E+ 


6.8484 X 10-3 


; Ft 

3. 1614 X 10-*^ 


1452.2 


72.26 


72 






2.2784 X 10-3 


5. 495 2 X 10 ° 


348. 31 


156. 24 


143 




1.906 4 X 10-3 


3.912 8 X 10-" 


412. 77 


180. 09 


171 


filH "2 T T 

°^Cu H 


A E^ 


1.702 2 X 10-2 


3. 379 6 X 10 


44. 41 


20. 01 


18 


Itt 

H2 


e Zj„ 


2.996 6 X 10-2 


1. 971 9 X 10 


58. 15 


14. 65 


14 


2()2tj 81 d 

^>J^Hg^ Br 


A (^^Ej 


5.209 5 X 10-3 


4. 833 5 X 10 


117. 44 


59.50 


54 


Hg^H 


A 2tt 

^ n3/2 


2.0248 X 10-2 


2. 0804 X 10 


56.61 


19. 81 


18 


X ^E+ 


5.9845 X 10-2 


2. 1268 X 10-3 


17.70 


6.20 


5 


(Hg^H) + 


^ ^E+ 


2.6564 X 10-2 


4.7212 X 10-3 


11.88 


7. 02 


6 




yliE+ 


1.886 5 X 10-2 


2.375 9 X 10-3 


16.75 


10.09 


8 






3.603 9 X 10-3 


2.119 9 X 10-5 


177. 34 


88.36 


80 


XiE+ 


2.855 5 X 10-3 


4. 171 1 X 10-6 


399. 79 


143.3 


134 






9.28469 X 10-3 


1.60515 X 10-4 


64.45 


32.74 


29 




A ^E+ 


4.175 81 X 10-2 


4.696 93 X 10-3 


11.91 


6. 107 


5 




(X ^E)'' 


5.276 71 X 10-3 


4.053 54 X 10-4 


40.55 


27. 08 


24 


'U2 


S ^B„ 


1.01746 X 10-2 


2.36197 X 10-4 


53.13 


28.06 


25 


X ^E^ 


7. 375 58 X 10-3 


1.650 10 X 10-"' 


200.98 


59.14 


50 



Table 2: continues 



molecule 


state, line 




\ye\ 


P 


3 


[n_] 


(24)MglH 




2.106 04 X 10-2 


1.002 87 X 10-^ 


81.53 


21.00 


20 


({24)Mgll^ + 




6.003 35 X 10-3 
1.78140 X 10-2 


3.178 25 X 10-4 
3.008 32 X 10-4 


45.80 
47.08 


29.33 
20.05 


26 
18 


p)Mg2H) + 




4.24725 X 10-3 

1.328 88 X 10-2 


1.43207 X 10-4 
1.36149 X 10-4 


68.23 

69.98 


43.46 

28.25 


38 

26 




a in. 


8.392 71 X 10^3 
9.51197 X 10-3 
7.559 65 X 10-3 


1.05646 X 10-3 
1.71189 X 10-^ 
2.06204 X 10-4 


25.12 
197. 34 
56.86 


16.62 
47. 80 
32.62 


14 
46 
29 






9.583 28 X 10-3 
7.310 65 X 10-3 


2.22122 X 10-4 
1.812 26 X 10-^ 


54.79 
191.80 


29.28 
58.91 


26 
56 


2^Na2 




5.09169 X 10-3 
4.55944 X 10-3 


7.561 19 X 10-^ 
1.69566 X 10-^ 


93.90 
198. 28 


52.31 
82.36 


47 
76 


O'J _ T 

^•^Na^H 


A 


7% 

1.74179 X 10-2 


1 

6.342 56 X 10-4 


32.42 


16.63 


15 


Na^^K 


A iE+ 


4.258 94 X 10-3 
1.092 05 X 10-3 


9.905 09 X 10-^ 
4.87163 X 10-^ 


82.04 
116.98 


50.14 
83.95 


44 
75 


14]N^16o 


^ 'ni/2 


6.106 36 X 10-3 
7.33842 X 10-3 
7.33707 X 10-3 


1.180 79 X 10-4 
6.303 58 X 10-^ 
6.30242 X lO-'^ 


75.14 
1028.40 

1028. 49 


42.35 
67.24 

67.25 


38 
67 

67 


^^02 




1.142 56 X 10-2 
9.736 95 X 10-3 


5.358 67 X 10-4 
7.503 38 X 10-6 


35.27 
298. 08 


20.45 
48.95 


18 
48 




. 0.-. ] 


2.98475 X 10-2 


2.03427 X 10-4 


57.25 


14.66 


14 






3.592 89 X 10-3 


6.829 57 X 10-6 


312.43 


113.18 


106 


'tTT 

rb'^^Br 


— '^T 

A (^S) 


n 

2.622 95 X 10-3 


1 

1.836 07 X 10-4 


60.26 


42.11 


37 


PbSe 


L> 


2.783 61 X 10-3 


2.100 84 X 10-^ 


178.14 


98.41 


88 


85Rb2 


C 


1. 843 15 X 10 ^ 

-L . U / ijO A lU 


±.^^yuz X lu 


136. 80 


90.04 
98 fin 


80 


28si32s 


E 

D In 


3.469 64 X 10-3 
4.64844 X 10-3 


8.153 66 X 10-^ 
8.789 06 X 10-5 


90.43 
87.09 


56.77 
51.15 


50 
45 




E 


6.314 77 X 10-3 


1.036 27 X 10-4 


80.21 


43.69 


39 


Su'-''0 


D iE+ 


5. 286 65 X 10-3 


2.317 20 X 10-^ 


53.6-1 


34.22 


30 


SnS 


E 


3.694 29 X 10-3 


4.06711 X lO-'' 


128.03 


71.77 


64 


SnSe 


E 


3.916 58 X 10-3 


8.138 35 X 10-6 


286.21 


103.74 


97 


SnTe 


I 
B 


5.44188 X 10-3 
6.643 51 X 10-3 


1.306 05 X 10-5 
5.64481 X 10-5 


225.93 
108.68 


76.64 
50.87 


72 
46 


Sr^^F 


B 


3.80446 X 10-3 


1.43179 X 10-5 


215.78 


94.81 


87 


Tl^^Br 


A {^litr 


4.75443 X 10-2 


2. 031 02 X 10-3 


18.12 


7.29 


6 


(ZniH) + 




2.03549 X 10-2 


1.043 84 X 10-4 


79.92 


21.53 


20 



14 



Table 3: Parameters of the model for negative molecular constants Xg and j/g 



molecule 


state 




IVel 


P 


f 


j 


l^-^CslR 




2. 794 12 X 10-^ 


1.71569 X 10-^ 


19.71 


10.02 


20.94 


(39)k2h 




1.133 22 X 10-2 


2.745 61 X 10-^ 


49.28 


25.65 


51.75 




A''^+ 


0.123487 


1.785 27 X 10-2 


6.111 


2.324 


7.275 




a''j:+ 


7.739 98 X 10-2 


6.520 91 X 10-^ 


10.11 


4. 084 


11.80 


^^Na^H 




1.74179 X 10-2 


6. 342 56 X 10-^ 


32.42 


16.63 


34.30 


Rb^H 


A'^+ 


1.676 21 X 10-^ 


6. 909 24 X 10-^ 


31.06 


16.441 


32.35 
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